use state.dta
frontier lnc_ lny1_ lny2_ lny3_ lnp1_ y1y1 y2y2 y3y3 y1y2 y1y3 y2y3 p1p1 y1p1 y2p1 y3p1 t tt y1t y2t y3t pt br r ry1 ry2 ry3 rp crr slr mac, distribution(tnormal) cm(r0 t0 tt0) cost nolog 
 predict e, te
 generate costeff = 1/e
